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Calculations relating to two experiments that demonstrate coherent control of preformed 
rubidium-85 molecules in a magneto-optical trap using ultrafast laser pulses are presented. The 
two experiments are a step towards the stabilization of ultracold rubidium dimers using ultrafast 
lasers. In the first experiment, it is shown that pre-associated molecules in an incoherent mixture 
of states can be made to oscillate coherently using a single ultrafast pulse. A novel mechanism 
that can transfer molecular population to more deeply bound vibrational levels is used in the sec- 
ond. Optimal parameters of the control pulse are presented for the application of the mechanism to 
molecules in a magneto-optical trap. The calculations make use of an experimental determination 
of the initial state of molecules photoassociated by the trapping lasers in the magneto-optical trap 
and use shaped pulses consistent with a standard ultrafast laser system. The experiment's purpose 
is to demonstrate and evaluate the use of ultrafast shaped pulses to manipulate ultracold rubidium 
dimers with a view to the eventual stabilization of the molecules. 

PACS numbers: 32.80.Qk, 34.50.Rk 



I. INTRODUCTION 

The success in cooling atomic gases to quantum de- 
generacy and the resulting rich phenomenology has 
prompted efforts to cool molecular gases to similarly cold 
temperatures. Two alternative approaches have been 
demonstrated that create cold molecules. The first ap- 
proach relics on methods that directly cool preformed 
molecules, but the lack of a closed cooling cycle pre- 
cludes the use of traditional atomic laser cooling tech- 
niques. As a result, molecules have not been directly 
cooled to the sub-microkelvin temperatures needed for 
quantum degeneracy. By contrast, techniques associat- 
ing precooled atoms retain the very cold ccntre-of-mass 
motion of their parent atoms, but the stabilization of the 
molecules in the lowest rovibrational levels is an ongo- 
ing challenge. To achieve a stable ultracold gas, the cold 
molecules must be transferred to their vibrational, rota- 
tional and electronic ground state. Spontaneous decay 
must be avoided if the very coldest temperatures are to 
be achieved since spontaneous emission heats the centre- 
of-mass motion of the molecules. 

Several recent efforts have been successful in stabilizing 
ultracold molecules. The majority of these rely on stim- 
ulated Raman adiabatic passage (STIRAP) [H, Typ- 
ically, highly vibrationally excited molecules are created 
using a magnetic field sweep across a Feshbach resonance 
which associates molecules into a highly vibrationally ex- 
cited state. STIRAP then transfers these to more deeply 
bound states using pairs of pulses. The process is very ef- 
ficient, and although there are constraints on the pairs of 
levels that population may be transferred between, mul- 
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tiple STIRAP transitions may be made that can populate 
the vibrational ground state. Notable demonstrations of 
this method include the stabilization of rubidium dimers 
into their triplet vibrational ground state with a single 
STIRAP pulse pair 0], the stabilization of KRb dimers 
into their singlet and triplet ground states Q , the trans- 
fer of caesium dimers using STIRAP to the seventy third 
vibrational level 0, then to the vibrational ground 
state 0- 

By varying a magnetic field and using a series of stim- 
ulated radio-frequency transitions, it is also possible to 
manipulate diatomic molecules but without exor- 
bitantly large magnetic fields, it is difficult to transfer 
molecules to any vibrational state bound by more than 
a few GHz. This technique can be used to achieve a pre- 
cise initial state for either STIRAP or some other control 
mechanism. 

A second route to the vibrational ground state of the 
caesium dimer has recently been demonstrated In 
this scheme, a continuous wave laser is used to photoas- 
sociate molecules into deeply bound vibrational states. 
Broadband pulses are used to excite the deeply bound 
caesium dimers. The pulses are shaped spectrally to 
avoid exciting the vibrational ground state. Over several 
cycles of pump and decay, population begins to accumu- 
late in the vibrational ground state. The disadvantage 
of this scheme is that it relies on spontaneous emission, 
and therefore the centre of mass motion of the molecules 
is heated. Also, the scheme relies on a fortuitous align- 
ment of excited state potentials in Cs2 for an initial injec- 
tion of population into deeply bound vibrational states. 
Such a mechanism has not been demonstrated in Rb2, 
and without some way to populate the lowest vibrational 
states, the pump-decay cycling results in the dissociation 
of the dimer with negligible accumulation in any vibra- 
tional state. 
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A complementary technique to STIRAP is control us- 
ing ultrafast laser pulses. Whereas STIRAP in general 
avoids spontaneous decay by avoiding populating decay- 
susceptible states, short pulse techniques transfer popula- 
tion on a timescale shorter than the decay time of the ex- 
cited molecule. The technique uses laser pulses with du- 
rations between tens of femtoseconds to several picosec- 
onds, and benefits from the relative case with which such 
pulses can be shaped spectrally. On paper, the achieve- 
ments of short pulse control on ultracold molecules have 
been modest. Photoassociation by ultrafast laser pulses 
has been demonstrated and the coherent control 

of photodissociation [illil of preassociated Rb2 has 
also been demonstrated. Molecular stabilization using 
ultrafast lasers has not yet been demonstrated. 

Short pulse control, unlike STIRAP, can address many 
transitions at once. This allows for a coherent final state 
to be chosen freely, allowing, for example, the creation of 
shaped wavcpackcts. If short pulses were to be used opti- 
mally and shaped with sufficient precision, they could be 
used with far more versatility than STIRAP. By contrast, 
STIRAP must be tailored to a specific transition, and 
each population transfer requires two phase-locked lasers. 
An experimental setup that achieves transfer across one 
transition using STIRAP is not able to address any other 
transition without modification. Short pulse control, by 
contrast, could change the final occupation of vibrational 
states in real time and with a far less restricted range of 
output states. Appl ying an adaptive technique such as 
a genetic algorithm [13| to the pulse shaper could then 
be used to create a desired superposition of states. Short 
pulse control, coupled with an adaptive system, could be 
used to manipulate molecules for which detailed spectro- 
scopic information is not present. For these reasons, it 
may be argued that ultrafast shaped pulses will have an 
important role in future developments in this field. 

A scheme has been proposed [13, [3] that stabilizes 
loosely bound *^Rb2 molecules by transferring them to an 
attractive excited state where they oscillate. A short time 
later they are transferred back to the electronic ground 
state. Optimal control pulses for such a 'pump-dump' 
scheme have been suggested [l6l |. 

Although it may be shown theoretically that a suf- 
ficiently tailored control pulse could transfer population 
to the vibrational ground state with high efficiency, there 
is uncertainty about how the calculations relate to ex- 
periments where the parameters may be controlled with 
finite precision. In this paper the gap is narrowed by the 
suggestion of two experiments and realistic calculations 
relating to them, that demonstrate several of these prin- 
ciples on a readily available source of molecules. The two 
experiments represent a step forward in the experimental 
application of control using ultrafast lasers to ultracold 
molecules. 

In order for this pump-dump scheme to work, the prin- 
ciples of its operation have to be shown to work: the cre- 
ation, detection, and manipulation of a suitable initial 
state; the creation, observation, and control of the os- 



cillating intermediate step; and the measurement of the 
final state. 

Unlike previous work, the calculations here relate 
to two experiments manipulating preformed molecules 
whose initial state was determined by analysis of ex- 
perimental data. The experiments are carried out on 
a dilute gas of ^^Rb atoms trapped in a magneto- 
optical trap (MOT). The gas cloud contains a reservoir 
of loosely bound ^^Rb2 molecules due to photoassocia- 
tion by the MOT trapping lasers and three-body recom- 
bination H, ■ The experiments manipulate these 
molecules. The parameters of the control pulses were 
chosen to correspond to a standard ultrafast laser sys- 
tem. 

The first experiment creates a coherently oscillating 
excited wavepacket. The experimental demonstration of 
this would open the door to wavepacket shaping and fo- 
cussmg [13, m, m, that is essential to the pump-dump 
scheme and other coherent control processes using ul- 
trafast lasers. A short laser pulse is used to excite the 
ground state (5S -I- 5S) ^^Rb2 molecules creating a co- 
herently oscillating wavepacket in the excited electronic 
states associated with the 5S -I- 5Pi/2 asymptote. The 
calculations presented here show that coherent oscilla- 
tions may be caused in the excited state despite the in- 
coherence of the initial state. 

The second experiment demonstrates the population 
transfer of loosely bound ground state dimers to more 
deeply bound vibrational levels and is a second way to 
demonstrate coherent control using ultrafast lasers. The 
process of optimization of the control pulses and the de- 
tection of the resulting molecules are similar to the equiv- 
alent processes in the pump-dump scheme. The mecha- 
nism used to increase the binding energy of the molecules 
may be applied to other molecular species. 

In Section ini the model used to predict the behaviour 
of an atom pair under the influence of a time dependent 
electric field is presented. The initial state of molecules 
photoassociated by the MOT trapping lasers is given. 
Sections mil and IIVI detail the experiments and give pre- 
dictions of their outcomes and, in the second experiment, 
the optimal pulse parameters. 



II. THEORY 

A. Model / Excited state dynamics 

A model containing eight electronic states is used to 
calculate the time dependent wavefunction of the dimer 
during and after its interaction with a time dependent 
laser field. The model Hamiltonian takes the form 



Hj.k = -S 



2n dm 



where j and k are the indices of the eight channels. Vj (R) 
is the potential energy of a given electronic state at a 
given intcrnuclcar separation, R. Dj ^ is the electronic 
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FIG. 1: (Color online) Eight electronic states are used in the 
majority of the calculations presented for reasons discussed in 
the Appendix. The eight electronic states in the model used in 
this paper have four distinct symmetries. Each symmetry has 
an associated potential energy curve. The potential energy 
curves for these four symmetries are shown here. They are 
taken from references [2l|, [H, [H, [l!] . 

dipole moment of the transition from state j to k. ^(t) 
is the time varying electric field of the pump pulse. Wj^k 
is the spin-orbit coupling between electronic states. 

The eight electronic states (see Appendix) consist of 
two ground states with (a)'^S]+ symmetry (with zero or 
one unit of spin angular momentum projected onto the 
internuclear axis.) and six excited states associated with 
the 5S + 5P asymptotes, of which three have ^Hg symme- 
try, two have '^E^ symmetry, and one has ^Ilg symmetry. 
These eight states form a set that are not coupled by ei- 
ther dipole transitions or spin-orbit coupling to any other 
states outside the set that are associated with the 5S -I- 
5S or 5S + 5P asymptotes. For reasons discussed in the 
Appendix, this eight state model was found to give the 
same results as the full 15 state model. 

Each symmetry has a distinct potential energy curve. 
The excited state potential energy curves for the dy- 
namics calculations are taken from ab initio calcula- 
tions [U, [23. Their long range behaviours are modi- 
fied to match published dispersion coefficients [l^. The 
ground state potential energy curves are taken from ref- 
erences [Ij] . Spin-orbit couplings and dipole transi- 
tion strengths are assumed to be constant and take their 
asymptotic values. These were inferred from the energy 
difference between the 5P3/2 and 5Pi/2 states and their 
lifetimes [26j. The potential energy curves are shown in 
Fig. [TJ The model neglects the rotation of the inter- 
nuclear axis which takes place on longer timescales than 
the experiments. The angle between the internuclear axis 
and the polarization of the laser field is fixed in each cal- 
culation. 

The initial state for the calculations is a density ma- 
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FIG. 2: (Color online) The relative occupation of triplet 
bound states in a *^Rb MOT's molecular reservoir as a func- 
tion of binding energy. The inset shows the potential en- 
ergy curves of the Og state (upper) and ground triplet state 
(lower). The MOT trapping lasers excite atom pairs at long 
range. Population is concentrated at the inner turning point 
of the Og potential where it can decay to loosely bound ground 
triplet molecules with the shown distribution. The relative 
occupation of each vibrational state was estimated by con- 
sidering its Frank-Condon overlap with a highly vibrationally 
excited Og state. The distribution is insensitive to the choice 
of vibrational state in 0^ as long as it is chosen close to dis- 
sociation. An RKR potential energy curve [2^ was used for 
the calculation of the state. The distributions agree with 
the experimentally measured distribution. 

trix, diagonal in the model Hamiltonian's eigenbasis. The 
diagonal elements are described in Section lll Bl The evo- 
lution of the density matrix in time was calculated us- 
ing a Chebyshev propagator on an analytic mapped grid 



B. Initial state 

Two mechanisms create molecules in a ^^Rb MOT: pho- 
toassociation by the trapping lasers, and three-body re- 
combination [l8|. At high trapping laser intensities and 
lower MOT densities, photoassociation dominates p^ . 
In the experimental setup described in this work, the oc- 
cupation of vibrational levels in the ground triplet elec- 
tronic state was found to match that predicted assuming 
photoassociation to be dominant. 

The trapping lasers in a MOT photoassociate unbound 
atom pairs creating a reservoir of loosely bound ground 
state molecules. The MOT consists of two lasers tuned 
to the D2 (5S — 5P3) transition in rubidium. To pho- 
toassociate, these must excite atom pairs to electronic 
states which are both attractive and associated with the 
5S -I- 5P3 threshold. Of the states that satisfy both cri- 
teria, photoassociation to the 0~ state is dominant. This 
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mechanism populates the top ten ground triplet vibra- 
tional states with a binding energy of up to 4.5 cm~^. 
The occupation of each vibrational state in the triplet 
ground electronic states was estimated by calculating 
their Frank-Condon overlap with a highly vibrationally 
excited 0^ state. 

Since the trapping lasers responsible for the photoas- 
sociation are tuned to within a few MHz of the dissoci- 
ation threshold, the 0~ vibrational state responsible for 
the photoassociation must be close to the dissociation 
threshold. It was found that the distribution of popula- 
tion in the ground state vibrational levels resulting from 
the decay of an 0~ state was insensitive to the choice 
of 0~ state, as long as it was close to the dissociation 
threshold. 

An RKR potential energy curve [1^ was used for the 
initial state calculation. The occupation of ground state 
vibrational states as a function of binding energy is shown 
in Fig. [2] with the relevant potential energy curves. The 
number of molecules in this distribution may be increased 
by using a continuous wave photoassociation laser tuned 
to photoassociate unbound atom pairs to the same 0^ 
state. 

The occupation of each vibrational state may be mea- 
sured by molecular spectroscopy experiments using a 
narrow-band laser tuned to around 14500 cm~^ l3, 29 1. 
In an accompanying experimental work to this paper ll| , 
a pulsed dye laser was scanned over a wavelength interval 
from 14200cm~^ to 14600cm~^. The number of molec- 
ular ions produced (via resonantly enhanced multipho- 
ton ionization through the (2)^I]+ state) as a function of 
wavelength was used to infer the population distribution 
of the vibrational levels of the molecular ground state. 
Good agreement was found between experiment and the- 
ory, suggesting that the detectable molecules were indeed 
created by photoassociation to the 0~ state. 
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FIG. 3: (Color online) The spectrum of the excitation pulse 
for the oscillation experiment in relation to the adiabatic po- 
tential energy curves, showing the two 5S -|- 5P asymptotes. 
The pulse is spectrally cut at 374 THz to prevent excitation 
of atoms to the 5Pi/2 state, and to prevent scattering pairs 
from being excited. 




III. COHERENT OSCILLATIONS 
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The reservoir of preformed molecules within a *^Rb 
MOT are an accessible basis for a coherent control ex- 
periment. However, they are initially in an incoherent 
mixture of the top ten vibrational states. It is not clear 
therefore whether or not they can be manipulated coher- 
ently. 

The model described in Section Hi Al was initiated with 
the density matrix described in Section lll Bl A time vary- 
ing electric field was used that approximates an exper- 
imentally realizable laser pulse: a Gaussian pulse with 
a centre frequency of 375 THz (799.2 nm), a full width 
at half maximum intensity of 8 THz (17 nm), and spec- 
trally cut to remove intensity at a higher frequency than 
374.7 THz (800 nm), removing spectral intensity at the 
5S — 5Pi/2 transition. The spectral cut prevents atoms 
from being excited in the MOT which depletes the MOT 
and floods the ion detectors with atomic ions. Having the 
spectral cut at 800 nm rather than 794.7 nm (the transi- 



FIG. 4: (Color online) The average internuclear separation 
and population density in the excited electronic states after 
being excited by a spectrally-cut short Gaussian pump pulse. 
Higher population density is represented with darker colors. 
Oscillations with a period of around 5 ps can be seen. The os- 
cillations persist with no appreciable decay to at least 250 ps. 
A fractional revival of the vibrational wavepacket is evident 
at 230 ac- 



tion wavelength) also prevents the photoassociation of 
unbound atom pairs to long range molecules which os- 
cillate more slowly. The spectrum of the cut pulse as 
well as the potential energy curves of the excited states 
are shown in Fig. [3l Calculations were performed at a 
range of fluences up to 300 Jm~^ with fluences in both 
the linear and non-linear regimes. 



5 



Results 



For all fluences, population is excited predominantly 
at an internuclear separation of 20 uq. It then oscillates 
with a period of around 5ps. The oscillations persist to 
the end of the simulation at 250 ps. The population den- 
sity and average internuclear separation with an incident 
pulse flucncc of 20 Jm~^ are shown in Fig. 21 Oscillations 
could be detected if a position sensitive measurement of 
the excited state population could be made. 



B. Detection of oscillations 
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The ionization cross-section has, in some molecular 
species, been shown to be dependent on the internuclear 
separation of the dimer [13, HH- It is therefore possi- 
ble that a suitable ionizing probe pulse can be used to 
observe the oscillations: The molecular ion signal would 
vary as a function of the time delay between the pulse 
causing the oscillation and the probe pulse, revealing de- 
tails of the oscillations. 

The experiment these calculations imply is similar to a 
previous experiment in which a magneto-optical trap 
was illuminated by a cut Gaussian laser pulse, and then 
probed by a time delayed probe pulse. The molecular ion 
signal observed consisted of molecules that had been pho- 
toassociated by the control pulse rather than preformed 
molecules. By contrast, the experiment here has a control 
pulse which is spectrally cut further to the red to extin- 
guish this second source of molecular ions. A continuous 
wave photoassociation laser tuned to a transition from 
unbound atoms to bound molecules in the 0^ electronic 
state associated with the 5S -|- 5P3/2 asymptote would 
increase the number of molecules in the initial state used 
here. This would also help to reduce the relative contri- 
bution to the molecular ion signal of molecules photoas- 
sociated from unbound atoms by the control pulse. 



FIG. 5; (Color online) The ground and excited states in the 
rotating wave approximation, showing how the mechanism de- 
scribed in Section [IV Al works. The ground state (solid line) 
is coupled by a strong electric field to an excited state (thick 
dashed line) . The coupled Hamiltonian can be diagonalized at 
each internuclear separation, resulting in two adiabatic poten- 
tial energy curves exhibiting an avoided crossing (thin dashed 
line). Near the crossing, both states are more attractive than 
the ground state. Population is not transferred to the excited 
state. This gives any population near the crossing a momen- 
tum kick inwards. 



Parameter 



Value 



Fluence before shaping 
Fluence after shaping 
Delay 

Blue cut (pulse 1) 
Red cut (pulse 1) 
Blue cut (pulse 2) 
Red cut (pulse 2) 
Chirp 1 
Chirp 2 



100 Jm-2 
6.17 Jm-2 
14.1 ps 

369.8 THz (810.5 nm) 
368.0 THz (814.4 nm) 
370.3 THz (809.4 nm) 
365.8 THz (819.3 nm) 





TABLE I: Pulse parameters for a pulse pair that achieves 
around 15% transfer from the ninth highest vibrational state 
to the twelfth highest in the ground triplet electronic state. 



IV. POPULATION TRANSFER 



Mechanism 



The reservoir of molecules may be expected to have 
predominantly symmetry. Pump-dump schemes as 
mentioned above are unlikely to work for this symmetry 
in rubidium due to the lack of any suitable excited po- 
tential energy curves. Instead, a scheme is proposed here 
that does not rely on an alignment of excited potential 
energy curves. The scheme increases the binding energy 
of the molecules using a dynamic Stark shift. Calcula- 
tions presented here suggest that the scheme is robust 
enough with respect to the control field that it is exper- 
imentally accessible. 



A mechanism is used to increase the binding energy 
of the dimer. It makes use of a momentum kick as dis- 
cussed in references [32l . Issj . The molecules are initially 
in a stationary wavefunction with an electronic configu- 
ration that has very weak attraction between the atoms. 
An oscillating electric field is applied on a timescale much 
shorter than the oscillation time of the internuclear sepa- 
ration, but longer than the oscillation time of the valence 
electrons. This adiabatically mixes the initial electronic 
configuration with another electronic configuration with 
stronger attraction between the atoms as shown in Fig.O 
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FIG. 6: (Color online) The temporal and spectral intensi- 
ties of a pulse pair that achieves 15% transfer from the ninth 
highest to the twelfth highest vibrational level in the triplet 
ground state. In a similar scheme, two identical pulses may 
be used with the parameters of the first pulse here. This sim- 
pler scheme results in 10% transfer instead. The parameters 
of this pulse are given in Table |l] 



The population which was in the ground state follows 
the dressed initial state adiabatically and remains in the 
ground state after the pulse. Because the dressed state 
is transiently attractive due to the position-dependent 
Stark shift, the dimer is given a momentum kick inwards. 

In the scheme presented here, the molecules then 
bounce against the inner turning point of the ground 
state potential. They are given a further momentum kick 
inwards at a later time by a second control pulse. This 
slows them down, trapping them in vibrational states 
with lower energy than their initial state. The tech- 
nique is applicable wherever two electronic states with 
different potential gradients are coupled by an electric 
field. Although the calculations presented here relate to 
molecules produced in a MOT, the scheme could also 
be extended to the loosely bound molecules produced by 
sweeping a magnetic field across a Feshbach resonance. 

In this application of the technique, the mechanism re- 
lies on the interaction between the ground state and 
the excited state with 0~ symmetry associated with the 
5S -|- 5Pi asymptote. This was determined by running 
the simulations with further reduced models containing 
different combinations of states (see Appendix). 



B. Pulse optimization 

The transfer of population between an initial state and 
a more deeply bound target state was optimized. The 
initial state in the optimization was chosen to be the 
ninth highest vibrational state — which is predicted to 
have the highest occupation as shown in Fig. [2] — rather 
than a density matrix, in order to reduce the computa- 
tion time of the calculations. The target state was chosen 
to be the twelfth highest vibrational state. The twelfth 
highest state was chosen because its binding energy is 
sufficient to distinguish it from any of the initially pop- 
ulated states, but is similar enough to the initial state 
that a detection scheme that had been demonstrated to 
work for the initial state would be likely to work still. 

Two laser pulses are used as the control field. The 
pulses mimic an experimental setup where a Gaus- 
sian laser pulse [with a central frequency of 375 THz 
(799.2 nm) and a full width at half maximum intensity of 
8 THz (17 nm)] is split into two equal pulses. Each pulse 
is separated into its Fourier components and is spectrally 
cut, removing intensity above a threshold wavelength and 
below a second threshold. Separate chirps were applied 
to each pulse. The pulse pair therefore has eight param- 
eters: the delay, the wavelengths of the two blue and two 
red cuts, the fiuence of the pulse pair before the spectral 
cuts, and the chirp of each pulse. The spectral cut has a 
sigmoid profile with a constant width in accordance with 
experimental capabilities. 

A genetic algorithm was employed to find a maximum 
in the transfer efficiency at a series of different fluences. 
These were examined and the most suitable — judged by 
having relatively high transfer efficiencies, low peak laser 
fields and low fluences — was taken as the optimal solu- 
tion. Low fluences allow the control fleld to be applied 
over a larger cross-sectional area of the gas cloud in the 
magneto-optical trap. High peak intensities are undesir- 
able because they increase the likelihood of multiphoton 
ionization during the control phase of the experiment, 
and also introduce Stark shifts that arc not taken into 
account in the models used here. 



C. Results 

At 100 Jm~^ (unshaped), 15% population transfer 
from the ninth to the twelfth highest vibrational state 
was achieved. The pulse parameters are listed in Table 
m and their spectra and temporal profiles are given in 
Fig.m 

The sensitivity of the transfer efficiency to the pulse pa- 
rameters was established by taking four two-dimensional 
cross sections through the eight-dimensional search space 
and plotting the transfer efficiency as a function of each 
pair of parameters as shown in Fig. [71 Although the inter- 
pretation of the plots is difficult, it can be seen that the 
scheme is most sensitive to the parameters of the first 
pulse and the delay between the pulses. The optimal 
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FIG. 7: (Color online) The efficiency of the pulse pair as a function of its eight parameters. The dependence on the fluence 
of the pulses and the delay between them is shown in the top left graph. The dependence on chirp of the two pulses is shown 
in the top right. The dependence on the centre frequency of the two pulses and their bandwidths is shown in the bottom two 
graphs. Negative bandwidth of a pulse results in the pulse being extinguished. Contours represent transfer efficiencies in steps 
of 2.5%. Higher efficiencies are represented by darker colors. Here, efficiency is defined by the transfer fraction from the 9th 
highest vibrational state to the 12th. 



delay corresponds to the time taken for the wavepacket 
to bounce against the inner turning point of the ground 
triplet potential and to return to the outer turning point 
of the target state, as can be seen from the population 
dynamics in Fig. [51 

The efficiency does not change dramatically as a result 
of nanometer shifts in the centre wavelengths of either 
pulse or in the bandwidth of the second pulse. The effi- 



ciency is sensitive to shifts in the bandwidth of the first 
pulse. It was found that the efficiency decreased roughly 
linearly with fluence for fluences less than the recom- 
mended fluence. The scheme is largely insensitive to the 
details of the second pulse, and it was found that replac- 
ing the second pulse with a replica of the first, still with 
14.1 ps delay, only decreased the efficiency of the scheme 
by one third. 
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FIG. 8: (Color online) The initial (dashed) and final (solid) 
populations of the highest vibrational states in the ground 
triplet before and after a control pulse pair. The initial pop- 
ulation results from a decay from a high lying 0^ state, and 
the final distribution shows how the pulse pair transfers pop- 
ulation to more deeply bound states. 



The bandwidth of the pulses used in the computation 
were limited to more than one nanometer by the width 
and shape of the cutting function. Therefore, it is possi- 
ble that with a sharper cutting function, narrower band- 
widths for the first pulse might be more efhcient. 

The calculations were extended to model how each 
of the initially populated vibrational levels respond to 
the control field. An initial density matrix as described 
in section III Bl was propagated to obtain a density ma- 
trix describing the molecular reservoir after the control 
pulses. As for the pure initial state, a significant fraction 
of the population was transferred to vibrational states 
more deeply bound than the initial distribution. The oc- 
cupation of vibrational states before and after the pulse 
can be seen in Fig. [51 The population density for the 
dimers as a function of time and internuclear separation 
is shown in Fig. [51 



D. Experimental considerations 
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This calculation implies an experiment in which a *^Rb 
MOT is illuminated by a shaped pulse pair and then 
probed via resonantly enhanced multiphoton ionization 
using a tuneable pulsed dye laser. A demonstrated de- 
tection scheme involves scanning the dye laser from 
14420 cm~^ to 14480 cm~^ to ionize the triplet ground 
state molecules. The number of molecular ions produced 
as a function of probe wavelength measures the occupa- 
tion of vibrational levels in the triplet ground state. 

The experiment depends on being able to make quite 
precisely shaped pulses and the fluences required are 
high, although experimentally possible. It is not clear 
how experimental considerations such as the presence of 
third order dispersion and other optical artifacts will af- 
fect the transfer efficiency. In turn, the calculations pre- 
sented here do not take into account losses to states as- 
sociated with the 5P -I- 5P asymptotes or multiphoton 
ionization which could reduce the number of molecules 
that are transferred to more deeply bound states. These 
losses are endemic to experiments using ultrafast laser 
pulses on molecules. Stark shifts caused by coupling to 
higher lying electronic states are also neglected. 



FIG. 9: (Color online) The population density as a function of 
time and internuclear separation. The ground (excited) state 
population is in the lower (upper) graph. Higher population 
densities are represented by darker colors. Population is in a 
stationary distribution initially. At Ops, it is given a momen- 
tum kick that results in a bounce against the inner turning 
point of the triplet potential at roughly 10 ps. At 14 ps it is 
stabilized by a second momentum kick. Although the wave- 
function in the ground state is strongly affected by the pulses, 
at no point is a significant population transferred to the ex- 
cited state. The parameters of the control field are given in 
Table m 



V. CONCLUSION 

Calculations relating to two experiments are presented 
here. The calculations show that the first experiment 
causes coherent oscillation in an incoherent mixture of 
pre-associated rubidium molecules with a typical period 
of 5 ps. 

A generalizable and experimentally realizable mecha- 
nism was found that increases the binding energy of the 
population of pre-associated molecules. A second exper- 
iment was proposed that uses this mechanism with opti- 
mized pulse parameters within experimentally realizable 
limits. 
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The two experiments would test and demonstrate spe- 
cific techniques such as shaped control pulses and time 
resolved detection of molecular oscillation in ultracold 
molecules, processes that will play a vital role in pump- 
dump schemes to stabilize ultracold molecules. More 
generally, these experiments are conservative ways to 
test the fundamental principle that shaped ultrafast laser 
pulses can be used to coherently manipulate ultracold 
molecules. 
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APPENDIX 

The rubidium dimer is modelled in this work by pro- 
jecting the electronic state onto a finite number of basis 
states but allowing the nuclear motion to evolve freely. 
The electronic basis states may be chosen to be eigen- 
states of the electron Hamiltonian with or without the 
spin-orbit contributions. The intricacies of Hund's cases, 
which specify the electronic basis set, are beyond the 
scope of this manuscript. The reader is referred to 
Ref. [H. 

The model used in the calculations presented here con- 
tains eight electronic states. A fuller model containing 
fifteen states and reduced models were also used for com- 
parison. 

Rb2 has four ground electronic states. One has sin- 
glet character and three have triplet character. Since the 
specific detection scheme used in this paper detects only 
triplet population, the calculations exclude the singlet 
ground state and the excited states not coupled to the 
ground triplet state. Of the triplet ground states ("^S^), 
one has 0~ symmetry, one has 1+, and one has 1^ sym- 
metry in Hund's case (c) notation, giving the symmetry 
under a reflection in a plane containing the two nuclei and 
the control field polarization axis as the superscripted 
symbol. All three triplet states may be populated by the 
Og decay mechanism. 

An electric field couples the triplet ground states to 
twelve excited states, given in Table HIl A selection rule 
for the interaction with a linearly polarized control field, 
and also for spin-orbit coupling is that the reflection sym- 
metry of the molecule does not change, and so molecules 
may be divided into those with negative and those with 
positive reflection symmetry. 

The dynamics in the excited states with 1+ and 1^ 
symmetry are identical as are the 2+ dynamics to the 2^ 



dynamics whereas the 0^ dynamics are not the same as 



Hund's case (c) 


Hund's case (a) 


Inclusion 
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no 
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no 
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gerade excited states 








no 


09" 




yes 




3y+ Irr 3tt 


no 


19- 


3y+ l]-r 3tt 


yes 


29+ 




no 


29" 




yes 



TABLE II: There are four ground electronic states of Rb2, 
and twelve excited gerade electronic states. The triplet states 
are coupled by an electric field to gerade excited states. If 
the field is linearly polarized, and the plane of refiection that 
defines the basis set is chosen to contain this field, then a 
selection rule exists that the reflection symmetry does not 
change. For reasons detailed in the text, the singlet ground 
state is excluded from the model as are all the states with 
positive refiection symmetry, giving the model eight electronic 
states. A fuller model containing all fifteen states (not the 
singlet ground state) was used for comparison. 



the 0^ dynamics. It was found empirically that omit- 
ting the states with positive symmetry did not change 
the excited state dynamics significantly, because the set 
of excited states with positive reflection symmetry do 
not contain any excited population on a potential energy 
curve that is not also found in the set of states with neg- 
ative reflection symmetry — both states with 0^ sym- 
metry are repulsive and so are only populated by laser 
intensity blue of the atomic 5S — 5Pi/2 transition which 
in both experiments was filtered out. It was also found 
that the mechanism used in the second experiment relies 
on the 0~ states and so only works for molecules with 
negative reflection symmetry. The states with positive 
reflection symmetry can therefore be omitted from the 
model because for both experiments the reduced model 
(8 states) gives the same dynamics as the full model (15 
states). 

Section flV Al discusses calculations making use of fur- 
ther reduced models. These contain the ground triplet 
states and the excited states with only one Hund's case 
(c) symmetry each. Whether or not each minimal model 
recovers the behaviour of the 8 state model determines 
which Hund's case (c) states are important in a given 
calculation. 
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